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Abstract: This paper investigates the influence of environmental noise on the 
characteristic timescale of the dynamics of density-dependent populations. Gen- 
eral results are obtained on the statistics of time spent in rarity (i.e. below 
a small threshold on population density) and time spent in commonness (i.e. 
above a large threshold). The nonlinear stochastic models under consideration 
form a class of Markov chains on the state space ]0, oo[ which are transient if 
the intrinsic growth rate is negative and recurrent if it is positive or null. In 
the recurrent case, we obtain a necessary and sufficient condition for positive 
recurrence and precise estimates for the distribution of times of rarity and com- 
monness. In the null recurrent, critical case that applies to ecologically neutral 
species, the distribution of rarity time is a universal power law with exponent 
—3/2. This has implications for our understanding of the long-term dynamics 
of some natural populations, and provides a rigorous basis for the statistical 
description of on-off intermittency known in physical sciences. 

Key words or phrases: Population dynamics - Stochastic non-linear differ- 
ence equations - Environmental stochasticity - Rarity - Ecological timescales - 
Power law - On-off intermittency - Markov chains - Martingales. 

1 Introduction 

Many of the traditional approaches in ecological theory are based on a paradigm 
that the ecological systems that we observe in nature correspond in some way 
to stable states of relatively simple ecological models. This is a view embodied 
in many classic approaches to ecological problems, ranging from the use of 
models to look at the coexistence of species to the interpretation of short-term 
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experiments [May [25] ]• This approach, however, has been questioned. As 
ecologists emphasize that long-term experiments yield different results compared 
with short-term experiments [Brown et al. |3J], theory must take the role of 
timescales into account and look beyond an emphasis on determining asymptotic 
stable states [Hastings |12|]. 

Of particular applied interest is the issue of identifying timescales that are 
relevant for population management. Study of pests or disease outbreaks often 
have the time spent by the system in quiescence as one of their focus. The con- 
servation of threatened species requires an understanding of how individual and 
environmental parameters affect the time spent by a population above some pre- 
defined critical threshold of abundance. Planning the exploitation of renewable 
resources in an uncertain world may also benefit from probabilistic estimates of 
time spent below and above predefined thresholds of ecological and economic 
interest. This points to a common mathematical problem: given information 
about the current population state, how can one characterize the distribution 
of time spent (or needed) by a population trajectory within (or to reach) a par- 
ticular domain of the population state space? Here we address this problem for 
a class of unstructured, density-dependent stochastic population models which 
includes the well-known Packer model. 

Theoretical ecologists have long used simple, unstructured models to inves- 
tigate general properties of population dynamics. These models, however, often 
miss essential features of natural populations, yielding predictions that turn out 
to be wrong even qualitatively [Durrett and Levin 1994 [3]]. The difference 
equation introduced by Ricker pi] makes a remarkable exception. Although it 
was originally introduced in the tradition of phenomenological modelling, the 
Ricker model has recently been re-derived from first principles accounting for the 
discreteness of individuals, the stochastic nature of birth and death events, and 
the spatial localization of interactions between individuals [Royama [M , Ripa 
and Lundberg |33| . Sumpter and Broomhead 36 ]. In a constant environment, 
the model reads 

X n +i = X n f(X n ) (1) 

where 

f{X n )=exp[r-aX n ], (2) 

and Xq > 0. X n measures the population size (or density at time n). The 
number r measures the per capita growth rate at low population density when 
density-dependent effects are negligible; thus, A = exp(r) is the population's 
geometric rate of increase from rarity. Unlimited growth is prevented by the 
term exp(— aX n ), where a > measures the intensity of negative interactions 
between individuals. This non- linear term typifies the mode of density de- 
pendence called overcompensatory, whereby xf(x) decreases toward zero as x 
becomes larger than some threshold (with undercompensatory density depen- 
dence, xf(x) would be monotonically increasing). 

Notation. Hereafter we set A = /(0) for the model (JJJ with general fucntion 
f(x) > 0; A is called the intrinsic growth rate. 
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The Ricker model shows a somewhat intricate transition from equilibrium 
and cycles to chaos as r increases [May |26|]. The natural timescale of the 
Ricker model in a constant environment is essentially determined by the intrinsic 
growth rate. When A < 1, the origin is the only stable equilibrium state of the 
system (hence lnX n goes to — oo and is thus transient), whereas when A > 1, 
there are two equilibria : the origin, always unstable; and a positive equilibrium 
which can be stable or not depending on the specific value of A (hence In X n 
is recurrent). Multi-dimensional Ricker models are also of interest to describe 
interacting populations of different species 

Randomness, however, is inherent to the real world. Perhaps the simplest 
way of incorporating stochasticity into (cf. [2HI) is to set 

X n+1 = X n f(X n )exp(Y n ) (3) 

where Y n , n = 0, 1, 2, • • • are i.i.d random variables with zero expectation; (X n ) 
is then a Markov chain (MC) on the state space ]0, oo[. This noisy Ricker model 
was first introduced in [Kornadt et al. 1991] and considered in the context of 
MC theory by [Gyllenberg and all |l()l l9"|]. The random variable Y n can be 
seen as an additive perturbation of the intrinsic growth rate, thus providing a 
model of environmental stochasticity jH] ; or as a perturbation of the population 
density which is 'felt' effectively by each individual in the population-a form of 
stochasticity that we call random heterogeneity. Random heterogeneity arises 
from differences in competitive abilities among individuals which may be rooted 
in intrinsic individual differences or in spatial variation in habitat quality |15|. 
In the first case a simple choice is to assume the Y n to be Gaussian whereas in the 
latter case, a typical choice for the law of Y n is the log-normal distribution . 
This model generalizes straightforwardly to the case of N species in interaction. 

Our goal in this paper is to study the time spent by the population within 
given ranges of density. Of particular interest are the time T e spent in a state 
of rarity defined by a small density threshold e > 0, 

T e := inf{n > 1 : X n > e} (4) 

when the current density Xq is smaller than e, and the time Tm spent in a state 
of commonness defined by a large density threshold M > 0, 

T M := inf{n > 1 : X n < M] (5) 

when the current density Xq is larger than M. More generally, we define the 
exit time of the Markov chain (X n ) starting at some Xq > from a domain A 
of population densities (in R + ) as 

T A = inf{n > 1; X n G A c } (6) 

so that T e = 7]o,e[ an d Tm = T]m,oo[- We further define the time of medium 
abundance as T[ e _ M ], and the time of escape from extremes as T^my- We then 
ask the following questions. 
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• What is the qualitative long-term pattern of population dynamics ? In 
mathematical terms, we ask whether the MC is recurrent, null recurrent or 
transient. This question was already adressed in [S] under the assumptions 
that the law of the random noise has exponential moments. Here we relax 
this assumption to consider moments of order strictly greater than one, 
thereby extending the scope of the model to random heterogeneity. 

• Then, we seek general results on the tail of the distribution of times of 
rarity or commonness which are independent of the current population 
density X . This means asking how the probabilities P(T e > n) (if X < 
e), P{T M > n) (if X Q > M) and P(T [e>M] c > n) (if X & [e, M}) behave 
when n goes to infinity. Almost equivalent is to ask for which p > 
is E[(T a )p] finite, given A =]0,e[, A =]M, oo[ or A =]Q, e[U]M, oo[ with 
X #A. 

• Another interesting feature of the distribution of the times of rarity or 
commonness concerns the dependence of their moments on the current 
density Xq, Therefore, we shall study the behavior of the pth moment of 
Ta, i-e. Ex (Ta) p for given A and p, as a function of the initial state Xq; 
namely, for given A =]0, e[ when Xq — > 0, or for given A =]M, oo[ when 

Xq — > 00. 

• How much do the preceeding results depend on the choice of the function / 
and the distribution of the noise ? 

• How do the results extend to multispecies interactions ? 

We shall consider the noisy dynamical system @ with a general continuous 
function / : K + ^]R + and independent copies (Yi,i > 0) of a variable Y such 
that E[Y] — 0. Throughout this article we shall assume the following : 

1. S[|y| 1 +' 5 ] < oo for some 5 > 0. 

2. The function f(x) > is continuous on [0, oo[ and there exist a > 0, r G K 
such that 

lim f(x)e ax - r = 1. (7) 

x — >oo 

The latter assumption says that at large density population growth is well 
approximated by the Ricker model. A crucial consequence is that the results 
will then only depend on whether the parameter A is strictly greater, smaller 
or equal to one, that is, whether the population is intrinsically growing, or 
declining, or neutral. For the Ricker model specified by @, these conditions 
translate into r < 0, r > and r = 0. 

We recall the meaning of basic terminology used for the classification MC 
processes. We say that the MC is transient if there is a positive probability 
that the time taken by the process, when initiated in any bounded segment 
A = [a,b] C]0, oo [ (with < a < b < oo), to return to A, is infinite; in other 
words, the probability that the MC never returns to A is positive. The MC 
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is recurrent if the return time to A is finite almost surely. The MC is positive 
recurrent if, in addition, the expectation of the return time to A is finite too 
and null recurrent if the expectation of this time is infinite. 

To help developing biological intuition for our analysis, it may be worth 
emphasizing that, heuristically, a (single) population's transience will indicate 
either extinction or escape toward infinite density. Finding that the distribution 
of time spent in a certain state has an exponential tail indicates that the popu- 
lation typically visits that states for short time only; in contrast, a 'polynomial' 
or 'heavy' tail tells that the population resides very long in that state. 

Recurrent notations in this paper will be P x for the law of (X n ,n > 0) 
starting from Xq = x, and E x for the expectation under P x . 

2 Time of commonness and time of medium abun- 
dance 

The distribution of the random time Tm of commonness (if M is fixed large 
enough and Xq > M) is described by the following Theorem l2.f I 

Theorem 2.1 (a) Asymptotics when X — > oo. For any X > M > we 

have: Ex q Tm < oo. Furthermore 

Ex T M -> 1, as X -> oo, (8) 
and Tm — * 1 o,.s. as Xq — > oo. 

(b) Asymptotics of Px (Tm > n) as n — •> oo, when Xo is fixed. 

From the fact that Ex Tm < oo, we have Px {Tm > n) — o(n _1 In -1 n) as 
n — > oo. 

Assume that in addition for some uq > E[e a ° Y ] < oo. Then there exists 
k < 1 such that for any Xq > M and any n > 1 

Px (T M >n)<C(X )K n (9) 

wra£/i some constant C(Xq). Moreover for all M large enough 

K < inf sup f(x) a E[e aY ] < 1. 

a>0 x>M 

If this injimum is reached for a — olq, then C(Xq) < (Xq / 'M) a ° . 

Part (a) of the theorem says that the negative density-dependence can pull the 
population out of extreme commonness in no longer than one step. Moreover, 
the distribution of the time of commonness is qualitatively little sensitive to the 
intrinsic growth A: its expectation Ex Tm is always finite. Consequently, the 
tail of the distribution of Tm, Px (Tm > n) decreases faster than (nlnn) -1 as 
n — ► oo for any current density Xo larger than M. In the case where the random 
noise has exponential moments (-E[e Qy ] < oo), this tail even decreases exponen- 
tially, as shown by equation 0) where the constant < K < 1 is independent 
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from the current density Xq, and the function C(X ) can be computed from the 
law of the environmental noise Y. Such estimate comes from the fact that the 
event of a large time of commonness commonness occurs only when an event of 
probability strictly less than one is repeated a number of times proportional to 
n (which is a standard large deviation estimate). 

The distribution of the random time of medium abundance Tu^] nas also 
an exponential tail. 

Theorem 2.2 Assume that Y is not uniformly bounded almost surely, i.e for 
all -oo < —M' < < M < oo, P(-M' < Y < M) < 1. Then, for all 
M > e > 0, there exists n(e, M) < 1 such that for all Xq £ [e, M], 

Px (T [e , M ]>n) < K(e,M) n . 

The hypothesis that the environmental noise Y is not uniformly bounded is 
satisfied by common choices of laws such as the normal or the log-normal laws, 
so that we did not try to weaken it. 

Let us now consider the two other random times T e and T[e,M] o > the behavior 
of which depends on whether A < 1, A > 1, and A = 1. 

3 Intrinsically declining populations 

When A < 1, the MC is transient with escape to 0. Thus, for any rarity threshold 
e > and any current density smaller than e, there is a positive probability that 
the time T t needed to cross the rarity threshold e is infinite IjlOl) . The behavior 
of the process is then trivial, as it follows from the law of large numbers [5]. 
The following formal statement is given only for the sake of completness. 

Theorem 3.1 Assume A < 1 and E[\Y\] < oo. 

Then the MC is transient. This means that for any compact subset A = 
[a, b] c]0, oo [ and any X £ A P Xo (T A c = oo) > 0. 

Furthermore for any e > and any Xq < e, 

P Xo {T e = oo) > 0. (10) 

4 Intrinsically growing populations 

When A > 1, the population MC is recurrent. Consider first the timescale of rar- 
ity: a small rarity threshold e > is fixed and the current density Xq is smaller 
than e. The expected time of rarity, Ex (T e ), is always finite. Furthermore, we 
show that for any X small enough, the expected time of rarity remains of the 
order of minus the logarithm of initial density, i.e. there are two positive con- 
stants c\ — ci(e) > and C2 = c 2 (e) > such that c\ < Ex (T e )/\ \nX \ < c 2 
for any initial current density Xq smaller than e. Similar to the case of the 
time of commonness (c.f Theorem 12. If) , the statistics of rarity time are given by 
Px (T e > n) going to zero faster than (n(lnn)) -1 when considering very long 
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time n, and this result holds irrespective of the current density Xq. Under some 
additional conditions on the environmental noise, namely i?[e^ ar ^] < oo for 
some a > 0, the distribution of time spent in rarity has an exponentially decay- 
ing tail. In other words it is controlled according to Px {T e > n) < C(Xo)p n 
for any time n > 1; here p < 1 is a constant independent of the current density 
Xo, whereas C(X ) depends upon Xq. 

Theorem 4.1 Assume A > 1. Then, the MC is recurrent, i.e for any bounded 
segment A = [a,b] c]0, oo[, Ta<= is almost surely finite: Px (Ta<= < oo) = 1 
for any Xq > 0. It is positive recurrent, i.e. Ex Ta<= < oo for any Xo > 0, 
iff E[e Y ] < oo. For any Xq < e we have Ex a T t < oo, so Px (T e > n) = 
o((n In n)^ 1 ) as n — > oo. We also have 

E Xo T e = 0(\lnX \), X ^0. (11) 

Moreover this order is precise in the sense that : 

lim inf > 0. (12) 

x ^o|lnX | 

Assume that for some a > Ee~ aY < oo. Then there exists p < 1 such that 
for any < X < e and any n > 1 

Px (T e >n) <C(X )p n (13) 

with some constant C(X ). Moreover, p < inf Q>0 e -" inf [o, £ ] fof(x) Ee -aY < 1 
If this infimum is reached for a = a^, then C(Xq) < {e/ Xo) a ° ■ 

Finally, we obtain information about the time Tj e ,M] c of escape from ex- 
tremes, that is the time needed for the population to reach a range of density 
[e, M] when initiated either below the rarity threshold e, or above the common- 
ness threshold M. In particular, we compute precise estimates for the latter's 
law, from which we derive the asymptotics of Px {Tum] c > n ) f° r l° n g times 
n. We show that this tail is of polynomial order ~ n~ a with degree given by 
the largest number a > such that -E[e" r ] is finite; thus, the tail depends on 
the law of Y via the value of a. This heavy tail arises from the very long time 
spent by the population regrowing from low density every time it falls short of 
passing below the rarity threshold e. 

Theorem 4.2 Assume A > 1. For any aM > r > e > 0, any Xq $l [e, M] and 

any a € M + , 

E[e aY ]<oo => ^[^.(lnr^c)* - 1 -^ ] < oo Vr? > (14) 
and conversely 

E[e aY ] = +*o E Xo [T£ M]c ] = +eo. (15) 
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In particular, if we let ocq = sup{a £ R + : E[e aY ) < oo}, then for any 77 > 
liminf n"°-" +1 P Xo (T [e . M] c = n) = 

n — >oo 

but 

limsupn Qo+? ' +1 P Xo (T [e!M] c = n) = +00. 

n — >oo 

This result calls for two further comments. First, compared to jf]], we show 
that the chain is positive recurrent if A > 1 under the sole condition that 
£[|y| 1+<5 ] < 00 (whereas assumed that E[e aY ] < 00 for some a > 0) and 
also specify the tail of Tu^e. Second, in lfP2|) . we in fact show (see (|57J) and 
H4UI) ) that there exists c, cf > such that for sufficiently small Xo 

c^supln/Cx))- 1 < Ex °f t] Y < cOn/CO))- 1 = c(lnA)- 1 . (16) 
x <e lne-lnAo 

This already indicates that the expectation of T e is a function of A and e which 
goes to infinity when A goes to one and e goes to zero, as stated in Theorem l5.ll 



5 Intrinsically neutral populations 

When A = 1, we assume that i?[|y| 2+(51 ] < 00 for some Si > 0. We also need 
to assume that f(x) —> 1 as x — * not too slowly: more precisely f(x) 
I + o(| ln:r|~ 1 ~' 52 ) as x — > with some small ^ > 0. This means that the 
population is neutral or almost neutral over a 'fair' range of small densities. 
If / were to go to one extremely slowly, one can intuitively guess that the 
effect of neutrality would be offset by the tendency of the population to grow 
or decline at very low density. Under these assumptions, the population MC 
is null recurrent. Thus, the expected time of rarity Ex T e is infinite, yet any 
initially rare population will escape from rarity in finite time with probability 
one: Px (T e < 00) = 1. We give a rigorous proof that the law of the time spent 
in rarity, i.e. Px (T e = n), is approximated by a power law with exponent —3/2 
as n —> 00. In contrast to the case A > 1, this law is universal, in the sense 
that it does not depend upon the distribution of the environmental noise. This 
behavior is the same as in the absence of population regulation (i.e. a = 0), 
when the population growth is driven by environmental noise only. 

Theorem 5.1 Let f(x) = l + o(| \nx\- 1 - 5 ' 2 ) as x -> and assume E[\Y\ 2+Sl ] < 
00 for some 61,62 > 0. 

Then the MC is null recurrent. This means that for any compact subset 
A = [a, b] c]0, 00 [ and any Xq $ A, Px (Ta<= < 00) = 1 but Ex Ta? = 00. 

For any p < 1/2 and for any 6 > there exists e(6) > such that for all 
X <e 

E Xo T e P = O(\\nX \ 2 P+ s ), X Q ^0. (17) 
For any p > 1/2, any e > and any Xq < e 

ExoT* = 00. (18) 
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In particular for any e > small enough, any Xq < e and any rj > there exists 
a constant C(Xq, e, n) > such that for all n > 1 

C(X ,e^)n-V 2 -< < P Xo (T £ = n) < C(X , e, r,) n -^ 2+r ' . (19) 

Compared to |S] we show that null recurrence holds also under assumption 
Ee aY = oo for all a > 0. All estimates given by the theorem are new. In 
particular we recover the predicted exponent —3/2 for Px„(T e = n), or equiv- 
alently 1/2 for Px (T e > n) (cf. 7 ). Thus, in such 'neutral' populations, the 
time spent in rarity is much longer than in commonness. The tail of this time 
is universal, i.e. it does not depend on the tail of the environmental noise Y as 
soon as Y has more than a second moment, and it is the same as for f{x) = 1 
provided the growth rate f{x) departs from one sufficiently slowly as density 
increases from zero. 

Ferriere and Cazelles |7] observed in their simulations that the exponent 
—3/2 arises for a suitable limit where A goes to one and e approaches zero. 
We here recover this result in the particular case where A is exactly one. The 
inequality (fH)|) also shows that if sup x<e \nf{x) goes to zero with e, the exponent 
has to be greater than 1. It is not yet clear how the —3/2 exponent could be 
recovered by taking such a limit. 

6 Two-species interaction 

The question of how these results extend to communities of interacting species 
opens a promising avenue for future research. As a first step in this direction, we 
report that the stochastic dynamics of two competing species with positive in- 
trinsic growth can be transient - a result contrasting with the monospecific case 
in which recurrence is guaranteed. Here, transience is understood in the sense 
that the two-species community will not visit all possible ranges of densities 
with positive probability. In fact, for some choices of the parameters, we show 
that there exists a parameter 8 > such that X^(X 2 )~ S is bounded below for 
all times with positive probability provided this inequality holds at time n = 0. 
However, the population X\ will undergo phases of commonness and rarity in 
finite time. Although the population density X 2 is somehow 'dominated' by 
X 1 , we could not resolve whether the population density X 2 was limited by 
some upper-boundary at all times with positive probability or not. 

A formal result can be stated by considering the two-dimensional Rickcr 
model with parameters 011,012,021,022 > 

Xl l+l = X l n e- a ^ xl ^ x ~ + ^ +Y ~ (20) 

^n+1 — 

(21) 

on ]0, oo[ 2 where (l^)„ e ^, (Y 2 ) n£ N are two independent sequences of i.i.d. ran- 
dom variables. 
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Theorem 6.1 Assume that E^Y 1 ) 2 + (Y 2 ) 2 ] < oo, n, r 2 > 0, n 7^ r 2 and one 
0/ </ie following conditions is fullfiled: 

1. ria 2 \ - r 2 an > 0, ria 22 - r 2 a\ 2 > 

2. ria 2 i - r 2 an > 0, ria 2 2 - r 2 a 12 < 

3. ria 2 i - r 2 an < 0, rio 22 - r 2 a i2 < 

Then the MC is transient. 

Furthermore, we can describe the following ways of escape to infinity. 

1. Assume that either condition (1) holds or r\a 2 \ — r 2 an > together with 
the assumption r 2 <r\. Let us fix e > such that (n — e)a 2 i — r 2 an > 0. Let 
M > 0. Let 

t m = inf {n : r 2 In X* - (n - e) In X 2 < M} (22) 
i/ r 2 In X X - (n - e) In X 2 > M. Then 

P(t m = 00) > 0. (23) 

Let Tj\(X l ) = inf{n > 0; X\ G ,4 C }, i = 1,2. T/ien, /or any L,e>0, 

Ex^niX 1 )} < 00, ExMX 1 ) I r M = 00] < 00, 

P Xo (T [eii]c (X 1 )Ar M <«)) = l. (24) 

2. Assume that either condition (3) holds or r\a 22 — r 2 a\ 2 < together with 
the assumption r 2 >r±. Let us fix e > such that (r\ + e)a 22 — r 2 a\ 2 < 0. Let 
M > 0. Let 

t m = inf {n : r 2 In X x n - (n + e) In X 2 > M} (25) 
r 2 In Xl - (n + e) In X 2 < M. T/ien 

P(r M = 00) > 0. (26) 

Let T A {X 2 ) = inf{n > 0; X 2 G A c }. T/ien, /or any L, e > 0, 

E Xa [T L (X 2 )} < 00, £* [T £ (X 2 ) I r M = 00] < 00, 

Px (T [(ii | C (I 2 )Ar M < oo) = 1. (27) 

We see that the two-population MC can be transient even for some r\ > 0, r 2 > 
while in the one-dimensional case it is always recurrent for r > 0. We conjecture 
that in the remaining case ria 2 i — r 2 an < and ria 22 — r 2 ai 2 > the chain is 
recurrent, as it can be suspected from the vector field of mean drifts, or from 
the isocline structure of the associated deterministic dynamical system. 
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7 Discussion 



The importance of identifying the biologically relevant timescales of ecological 
models, rather than focusing on their asymptotics, albeit emphasized only re- 
cently [Hastings has long been recognized in some areas of population 
biology. In the study of epidemic dynamics, for example, the relevant timescale 
can be the scale of a single outbreak, or longer if the question of interest is 
the timing between epidemics [Finkenstadt and Grenfell 8 ] . Another example 
arises in the study of plankton community dynamics in temperate regions, where 
seasonality can reduce the relevant timescale to less than one year, making the 
asymptotic behavior of models meaningless [Huisman and Weissing, ^H], |19|]. 

Environmental stochasticity can play a crucial role in determining the rele- 
vant timescales of a population's dynamics. For example, adding stochasticity 
to a metapopulation models could turn convergence to a stable focus equilib- 
rium into a kind of dynamics that would likely be identified as a noisy cycle 
by empiricists [Gurney et al. |28|]. implying that the deterministic equilibrium 
provides limited biological insight; instead, the period of the apparent cycle 
unravels a relevant timescale to interpret fluctuations in such a system. Even 
small noise can induce rapid fluctuations of large amplitude in population mod- 
els that would otherwise converge to stable equilibria [Higgins et al. H3|], e.g. 
by stabilizing population trajectories around chaotic repellors [Rand and Wilson 

ED]. 

Gyllenberg et al. JT] were the first to establish conditions under which 
recurrence and null recurrence occur in the Ricker model studied in this pa- 
per; they also addressed the case where the noise affects our parameter a, which 
represents a random perturbation of the carrying capacity of the habitat. Fager- 
holm and Hognas [Hj extended their results by considering the case where our 
parameter a is also affected by i.i.d. random perturbations. Our analysis offers 
a natural continuation of these previous studies that (i) assumes less stringent 
conditions on the moments of the random variables, (ii) provides precise and 
rigorous estimates of the laws of several characteristic times, (iii) paves the way 
for tackling similar issues in multispecific communities. 

Extension (i) is interesting as it broadens the scope of the model to the case 
where noise reflects random heterogeneity. Extension (ii) has several important 
implications. Building on previous heuristics and numerical studies, it yields a 
rigorous basis for the phenomenon of on-off intermittency. The notion of on- 
off intermittency refers to a certain class of burst-and-crash dynamics; it was 
first described in physical systems (Heagy et al. an d later found to be 

relevant in the ecological context (Ferriere and Cazelles [7|). We have shown 
(cf. Theorem 12. 2[) that for unbounded noise the distribution of time spent in a 
medium abundance has exponentially decaying tail. Therefore, the population 
spends most of the time fluctuating in states of commonness or rarity. The 
time spent in these states can be influenced by the tail of the environmental 
noise when A > 1. If A > 1, the distribution of the exit time of these extreme 
states can have heavy tail (cf Theorem l4.2f> even though the distributions of the 
time spent in either commonness or rarity have exponential time. If A — 1, the 
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population spends most of the time in rarity with occasional short outbursts 
(cf. Theorem I5.H . For such neutral species, the distribution of rarity times 
can be approximated by a power law with exponent —3/2; this approximation 
is universal in the sense that it does not depend upon the law of the noise, 
provided the latter has finite second moment. Ferriere and Cazelles 7 found a 
good fit of this approximation to real data on fish population dynamics. 

Many previous mathematical studies of population models based on stochas- 
tic difference equations trace back to the seminal work of Ellncr 5 . The first 
detailed study of a noisy version of the Ricker model similar to ours was carried 
out by Schaffer et al. [33] • Kornadt et al. [231 have considered a different model 
of environmental stochasticity and investigated how stability of the non-trivial 
equilibrium and period-two cycles of the deterministic Ricker map was affected. 
Similar questions were raised by Sun and Yang 37] for a model of noise reflecting 
random immigration. Hognas |17| studied the quasi-stationary distribution of 
the MC on a countable state space (branching process) that describes a popula- 
tion regulated by the Ricker mechanism with A > 1 . In Athreya and Dai |2] , the 
state space is continuous but the regulatory mechanism is given by the simple 
logistic map; they study convergence and the existence of a stationary measure. 
Ramanan and Zeitouni considered small noise operating additively on an 
iterative map of a compact interval, like the logistic map. In their case, extinc- 
tion (i.e. transience in this study) was possible even with A > 1, and a tension 
arose between this effect of stochasticity and the deterministic dynamics; for 
small noise, this tension resulted into large extinction times. There have been a 
number of large deviation analysis conducted on discrete-time models relevant 
for ecological applications. Morrow and Sawyer [27J studied the large deviation 
tail of the exit time from a neighborhood of a stable fixed point of a MC in the 
limit where the noise vanishes as time goes to infinity; Kifer |211 122"] studied 
in a broad generality the large deviation estimates for small perturbations of 
dynamical systems. 

Whether the universal power law found for a single neutral population exists 
for multi-specific communities warrants further investigation. From an ecolog- 
ical viewpoint, it would be interesting to elucidate the potential relationship 
between the exponent of such putative laws, and the structure of the network 
of interactions — how many species involved, who interact with whom, and how? 
The goal then might be to find some signature of a community's structure in 
the time series of those neutral or quasi neutral species (i.e. having a A close 
to one) exhibiting on-off dynamics. Another desirable extension should aim at 
incorporating temporal autocorrelation in the noise. Temporal autocorrelation 
in random environments can have dramatic consequences on the dynamics of 
density-dependent populations (Kaitala et al. [20]) Ripa and Lundberg 
Yet, the effect of autocorrelation at the characteristic timescales of rarity and 
commonness remains completely unknown. 
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8 Appendix : proofs of the theorems 

The proofs of all our statements rely on martingale techniques. 

Notation. Everywhere below we denote by Ex a x ; the conditional expectation 

over the sigma- field Ti — J-{Xq, . . . , JQ). 

8.1 Proof of Theorem IZTl 

First, we need to get an a priori estimate for Ex Tu (EEJ>- For any X large 
enough, define 

n 

lnX n = \nX +J2 ((M - InX^x) V (-1)). 

i=l 

Due to the assumption J7J on the function f(x), we may assume below that 
f(x) < e - ax+r+£ for all x > M and some e > 0. Then 

E[\nX n -\nX n ^ \ X n ^ > x] 

< E[(-aX n ^ +r + Y n + e)V (-1) | X n ^ > x] 

< (-ax + r + e + 1)P(Y > -1 + aX n - X - r \ X n - X > x) - 1 

+ E{Yl{ Y >-l+aX rl - 1 -r-e} I > x] 

< -l + E[Yl { Y>-l+ax-r-e}] = - 1 + o(l), X -> OO 

where we used the fact E[Y ly>-i+ a i-r-e] - * as x — > oo due to the assumption 
E\Y\ < oo. Thus for all x large enough 

E[\nX n - lnXn-! \ X n _! > x] < -1/2 

for all n > 1. It follows that 

£J[lnX nA T M | F n -x] < lnX„_i A T M - ~l{T M >n-l}, 

hence 

£x [ln^nA2v] < Sx [ln^ n -iAT M ] - 2 P x ( T M >n-l), 
and therefore 

1 " 

£x [lnX nAT J < -g P x (T M > i - 1) + lnX . 

But lnX„ > lnX„ and lnX n > lnX„_i — 1 > lnX„_i — 1 so that lnX nA T M > 
In M — 1 by the definition of Tm and In X n . Then 

1 " 

InM - 1 < --^2Px {T m > i - 1) +lnX 
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for all n > 1, whence 

oo 

E Xa T M = Px (T M > i - 1) < 2(lnX - In M + 1). (28) 

i=l 

Note that this inequality was obtained under the mere assumption that f(x) < 
e -ax+r+e f 0I ^ x > M and some e > 0. The next result requires the more 
precise estimate that for x > M with M fixed large enough 

e -ax+r-t < < e -ax+r+e 

for some e 6 (0, ir). We have the following renewal equation: 

Ex T M = Px„ (T M = 1) + E Xa [l Xl >M E Xl [T M + 1]]. (29) 

The first term in this sum is the probability that the set ]0, M] is reached in one 
step. Since f(x) < e - ax + r + e for x > M, we find that 

Px (T M = 1) = Px (*i < M) > P Xo (X e- aX " +r+c+Y < M) 

= P(Y <lnM -lnX + aX Q -r-e). (30) 

By the assumption i?|F| 1+5 < oo we have P(Y > t) < Ct^ 1 ^ 5 with some C > 
and alH > 0, which yields when In M — InXo + aXo — r — e > 

Px (T M = 1) > 1 - C(lnM - lnX + aX - r - e)- 1 - 5 . (31) 

Thus 

Px (T M = 1)->1, Xq^oo. (32) 
The remaining term is bounded according to 1|28[) by 

E Xo [l Xl>M E Xl [T M + 1]] < E Xo [l Xl>M {l + 2(lnX! - InM + 1))] 

< E Y {lY>inM-inX +aX - r -e(l+'2{Y+lnX -aX +r+e)-2\nM+2)). 

Again the assumption _E|Y| 1+<5 < oo and Chebychev's inequality imply, with 
A = hi M — lnX + aX Q - r - e, that 

E[l Y >lnM-lnXo+aX -r-e{l + 2{Y + lllX - aX + r + e) - 2 InM + 2)] 

= E[1 Y > A (3 + 2(Y - A))} < A- 5 E[\Y\ S (3 + 2\Y\)} 
< C'QnM -lnX a + aX a -r-e)- 5 

with some constant C > for all X$ large enough. Hence 

ExoIIx^mExATm + 1}} ->0, Xq^oo. (33) 

Combining J5?5J|, (J23) and we deduce that E Xo T M — ► 1 as X — > oo. Fur- 
thermore by (|3*TJl 

Px (2m > 1) ^C^ 1 - 5 
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with some constant C\ > 0. Then Tm — > 1 a.s. as Xq — + oo by Borel-Cantelli 
lemma. 

Let us now prove the estimate |(j5J. By Chebyshev's inequality for all a > 

Pxo(T M > n) = Px (ln^i > hiM Vt = 1, .. . , n) 

= P Xo (e°( ln x "- ln x °) > e a < ln M - ln x °\ X 4 > M Mi = 1, . . . , n - 1) 
S e u; ii[e ™ u 'J-{Xj>MVi=i,...,n-i}J 

(v IA/T\ar? ^a(ln Xi — In Xo) i tt- „a(ln — In Xi ) i 

- (A /A4) iijsfo 6 v ' *-{X 1 >M}&Xo,Xi e y 'L{X 2 >M} 



• • • 1 



{X n - 2 >M}E Xo ,...,X n 



,a(ln X„_ 1 -lnJC„_ 2 )l 

3 ; 1{X„_!>M} 



^X ,...,X„_i 

Note that for any i = 0, 1 

L{JCi>M} 



,a(ln X n -lnX„_i) 



,n — 1 



(34) 



*-{x i >M}&x ,...,x i e + < sup }(x) tie 

x>M 



Applying this inequality subsequently for i — n — l,n — 2, . . . , to the right- 
hand side of © we bound its left-hand side by (X /M) a K(a) n where k{q) = 
sup x > M f(x) a Ee aY . 

Our final task is to optimize this bound over a > 0. The function n(a) 
exists at least for a G [0, ao] as Ee a ° Y < oo. Furthermore k(Q) — 1 and 
k'(0) = sup x>M ln/(x) < 0, from where K,(a) < 1 for a > small enough. 
Thus k = inf Q> o n(a) < 1 and @ follows. 

8.2 Proof of Theorem 1231 

The proof of this theorem is straightforward since 

Px (T [6tM ] > n) = P Xo (ni< fc < n - 2 {X fc /(X fe )e yfe e [e,M],X k e [e,M]}) 

< Px„ (ni< fc < n - 2 {^ mi n(e, M) + e < Y k < X max (e, M) + AT}) 
= K (e,M) n - 2 



with K min (e, M) = -logm&x xe[ ^ M] (xf(x)), K max (e,M) = -logmm xe [ eyM ](xf(x)), 
and 

«(e, M) = P Xo (X min (e, M) + e<Y < K max (e, M) + M) < 1 

since we assumed that Y is not compactly supported and < min x6 r ej M] ( x f( x )) < 
max a . e[CjM ](a;/(a;)) < oo. 
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8.3 Proof of Theorems H~T1 and ICT 

Let us start by the proofs of and (|12f> . Assume that e is fixed small enough. 
We fix e > e and introduce X n = X„ A eo- Then 

£(mX n -lnl„_i | X„_!=a;) 

= £[(-m/(x) +y)l{Y<-ln( a; /( K ))+ln £ o}] ( 35 ) 

+ £[(lne -lna:)l { y > _i n(a/(x))+lneo} ]=ln/(0) + o(l)>]ii/(0)/2 ) x -» 0. 
Here we used the facts that ln/(0) > and, due to i£Y = and i?|F| 1+ ' 5 < oo, 

E[Yl{Y<-\n(xf(x))+\nc }} = E[Yl{ Y >- ln(x/(x))+ln e }] ~ * ^ ^ 0, 

and by Chebyshev's inequality 

E\Y\ 1+S 

P(Y > -ln(,/(,)) + lne ) < ( _ - ^ ^ = o((-ln,)-), 

as x — > 0. Thus lnX nA T" e is a negative submartingale and for e sufficiently small 

£x [lnX nATe ] > £ Xo [lnAViAT c ] + ^^Px (T e > n - 1) 
resulting, for all n > 0, with 

In /Yen n— ^ 

In e > E Xo [In X nATt ] > In X + ]T P A - (T £ > p) . (36) 

Hence, T e < oo almost surely as a consequence of Borel Cantelli's lemma. More- 
over 

l ^^-E Xo [T e ] <lneo-mX (37) 

which proves 

To get the lower bound on E Xo [T e ] claimed in (|12fl . first of all, starting from 
the rough estimate 

E[laX n - lnX n _i | X n -i <e]< E Y [{d + Y)]=d (38) 

with d = sup xe r e i ln/(x) > and proceeding in the same way as for (|36|l we 
obtain the upper bound 



E Xn \nX nATt < dY,Px (T e > i - 1) + lnX . (39) 

i=l 

Since T e < oo a.s. by (|1 lf> . then lnX„ A T e — > ln-Xr,, a.s. Note also that 

oo 

In X nATc | < = | In X k AT e - In X k -iAT e \ + In Xq. 



k=l 
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We have 

E(\lnX kATc -lnX fe _i ATe | | JF fe _ x ) 
= E(\lnX kAn -lnJffc_i AT .|lT.>fc-i l^fc-i) < (d + £|y|)l T .>fc-i. 
Then £J| lnX feATs - lnX fe _i ATe | < CP(T e > k - 1) with some constant C > 0. 
We have already shown that Ex T e = J2ilo ^ x o > «) < oo, then is finite. 
Hence, the dominated convergence theorem applies to the sequence In X n ^T t ■ 

Ex In X nA ,x t — > -Exo In Xr e • 

But from the definition of T e we have Ex Q ln-X^ > lne. Thus taking the limit 
as n — > oo in (|39[) we get the lower bound for Ex T € valid for all Xo small 
enough: 

{\ne-\nX Q )/d<E Xo T t (40) 

showing l|12|). 

The proof of H13[) is analogous to that of J3J. Namely, 

Px (T e >n) =Px (- In X t > -lne Vi = l,...,n) 

= P Xo (e a (- lnX " +lnX °) > e Q( ^ lne+lnXo) ,X < e, V» = 1, . . . , n - 1) 
< (e/X rE[e^' lnX ^ x "h { x t<i ^...,n-i}} (41) 

for all a > 0. We rewrite the right-hand side of (|41|l as a sequence of conditional 
expectations like in (|34fl . Each of them is bounded from above by: 

< l {Xt< e}f(X l )- a Ee- hY < e - Qinf i«-i ln fE[e- aY ] = p(a). 

Then the right-hand side of l|41l) obeys the upper bound (e/X ) a p(a) n with 
p(a) finite for all a > small enough, p(0) = 1 and p'(0) = — infro, e ] In / < 
for e small enough. Hence inf ct> o p(a) < 1 and (| 1 31) is proved. 

We are now ready to prove the recurrence of the MC. We denote by A the 
compact A = [e,M], then T[ ei M]° = inf{n > 1 : X n G A}. We fix e > small 
enough and M large enough such that 

P(Y>lnM- sup In f(x) -lne) =7 < 1. 

ase[Q,e] 

We prove that for any current density Xq G A c the time to reach [e, M] satisfies 

T[ etM ]c < 00 a.s. (42) 

Let T 1<M = inf{ra > : X n < M}, T lj( = M{n > T 1<M : X n > e}, . . . , 
T k ,M = inf{n > T fe _ M : X n < M}, T k . e = inf{n > T k . M ■ X n > e}. Then by 
© and (JTTJ 

Pxo(nr = i{Tfe,M < oo, T Ke < co}) = 1. (43) 

But 

^x (r[e,M]c > T*, e ) < P Xo (X Tl , c > M, X T2 t > M, .. . , Xr fc(e > M) 
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— E Xo [l{X Tl e >M,X T2ie >M,...,X Tk _ l t >M}Px Tk _ l f ( X T k , c > M)] 

< E x „(l{x Tl e >M,x T2 t >M,...,x T >M})P(y >lnM - sup \nf(x)-\ne) 

xE[0,e] 

<■< [P(Y > hxM - sup ln/(x) - lne)] fe = 7 fe -> 0, k -> oo, 

i£[0,e] 

where < 7 < 1. It follows that 

Px ( U^ = i (T [e , M] c < T fe , £ )) - lim Px (T [e , M ]c < T fc , e ) - lim 1 - 7 fe = 1, 

implying together with l|43|l that Tuj^-ia < 00 a.s. 

Let us show that B[e y ] = 00 implies Ex [TLjiil ] = 00 which proves that 
the recurrence of the MC is null. First of all, we remark that Ex 1h.Xt m = 
£;xo(ln^il{T flf =i})+- B ^o( m ^T JV /l{T M >2})- Here, due to the inequality e -<™+r-t> < 
f(x) < e ~ ax + r + 5 for x > M and some S > 0, we have: 

ExoilnXxl^^y] > E Xo [0-nX o - aX Q + r + Yj_ - 5)l{ Yl <lnM-r+aX -lnX n -8}] 

> -00 (44) 

but 

E Xo (\nX TM l {TM > 2} ) < E Xo (lnX 2 lT M >2) 

= E Xo [(In X - aX + 2r + 26 + Y 1 

n Y —aXo+T—5+Yi , v \-r , 1 
— aVS.0 e + r 2j- L {Yi>lnM-r+oXo-lnX -<5}J 

= -OO, (45) 



whence Ex [hiXT M } — —00. As we have £ , x [ m ^T M l{x TM >e}] > —00, then 
E Xo [lnX TM l { x TM<e }} =-00, that is 

E Xa [| In X Tm \1{x Tm <e}} = +00. (46) 

Let us also observe that by |Q there exists a constant C > such that for all 
X <e 

E Xo T e >C\\nX Q \. (47) 
We may now conclude from (|46|l and 147(1 that 

E Xo [T[e,M]°} > £ ! Xo-E ! X TM [r e l { x TM <e}] > #X [C| In X T M |1{X Tm <e}] = +00- 

(48) 

The positive recurrence of the MC under the assumption Ee Y < 00 has been 
proved in UJ. The proof in the genereal case of f(x) with ln/(0) > is com- 
pletely analogous: it uses the same test function g(x) — xli x > xo y+a \nx\l^ x<1 
(where xq is such that cixq = | lnxo|) for which the expansions (|50|l . (|51|l with 
a = l are valid. Therefore we do not repeat it. 
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Finally we turn to the refinements claimed in (|14f> and (|15|) . Let us first 
prove l|15(l . Observe that for any finite Xq ^ 0, -E[e Qy ] = +00 implies that 

£[l{x 2 <i}(ln^r] =+00. (49) 

Indeed, 

lnX 2 < {lnX -(aXo-r) + 26 + Y Q + r + Y 1 )-{aX e- (aXo - r ^- 5 )e Yo :=Z 1 -Z 2 

with < 00 since Yx,Yo have finite first moment, but E[(Z2) a ] = +00 

by our assumption. (|49|l follows readily. We now prove that there exists a a > 
and b a finite such that for any Xq < e 

E Xo [T?]>a a l {Xo < e} (ln-^) a + b a 

This is enough to get the desired estimate by starting at time n = 2 and trans- 
lating T c by 2 as in the proof of the null recurrence. Set S n = \nX n so that 

nAT e nAT, 

SnAT e =S +J2 ln /(^c) + Yk - 
k=0 fc=0 

We define M n — J2k=o^ k - Since d = sup^g^ ln/(a;) > 0, then ln/(Xfe)'s 
are non-negative in this sum and bounded from above by d. Then for any time 
n > 

S nA T e < S Q + d(n A T e ) + M„ A T e - 
In particular, since on T e < n, S nA T t > lne we get 

lne< (S + d(T e An)+M T , An ) 

so that for all neN 

l{T«<n} x [(lne - So) V 0] < 1 { t c <«} x [(d(T t An) + M Tc An) V 0]. 

Integrating the power a > of this inequality and using that for any a, b G K, 
any a > 0, \a + 6|" < 2 Q (|a| Q + \b\ a ), we obtain 

E[[(ln V 0] Q l m <„ } ] < £[|d(T e A n) + Mt £ a„)| q l { r e <n}] 

< 2 a d a E[(T e A n) Q ] + 2 Q £[|M TeA „n. 

By Burkholder-Davis-Gundy inequality (sec c.g 32 , Theorem 4.1 p. 160), since 
(M n )„>o is a martingale with < M > n — E[Y 2 ]n, we know that for any a > 0, 
there exists a finite constant C a such that 

£[|M T5A „| Q ] < C a E[{T e A n)f ] < C a (l + E[(T e A n) Q ]) 
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so that we get 

[(In -L.) V 0] a P(T t < n) - 2 a C a < 2 a (d a + C a )E[{T e A n) a ]. 

We can now let n going to infinity and use monotone convergence theorem with 
the fact that T e is almost surely finite to conclude 

[(In 4-) V 0} a - 2 a C a < 2 a (d a + C a )E[T*] 
Xq 

which finishes the proof of (|15|) with l|49|l . 

Finally it remains to prove (|14[1 . Let us introduce a positive function g(x) — 
\a~ 1 lnx\ a for x < xq and g(x) = x a for x > xq with xq chosen so that g is 
continuous. Here, a is such that lim^^oo e ax ~ r f(x) = 1, hence a~ 1 x~ 1 ln/(cc) = 
— 1 + 0(x~ 1 ) as x — > oo. Then by the assumptions = 0, -El^l < oo and 
Ee aY < oo we have the following asymptotic expansions 

E[g{X n+x )-g{X n )\X n =x] 

= E[(\a,- 1 \n{xf{x)) + a~ X Y\ a - z a )l{y<-in/(x)-ln(x/x )}] 
+ E[(x a f(x) a e aY - x a )l {Y ^ lnf{x) _ lD{x/xo)} } 

= x a (-a) — (l + o(l)) = -aa- 1 x a - 1 lTix(l + o(l)) 1 x -> oo (50) 
ax 

and 

E[g(X n+1 ) - g(X n )\X n = x] 

= ElQa-Hnixfix^+a-'Yl - ^ lnx\ a )l {Y< ^ ln f{x) . ln{x/xo)} } 

+ E[(x a f(x) a e aY - |a- 1 lnx| a )l {r >_ ln/(x) _ ln(K/xo)} ] 
= la" 1 lnx\ a E[(\l + (ln/Oz) + ^(lnx)- 1 !" - 1)1 {y< _ 1b/(i) _ 1dWio)} ] 

+ o(| ln^l"- 1 ) 
= |a- 1 lnx| Q a(ln/(a:))(lna;)- 1 (l + o(l)) 

= -aa- Q (ln/(0))|lna;| Q - 1 (l + o(l)), x 0. (51) 

Let us now start from some Xo ^ [e, A/] with e chosen small enough and M 
large enough. Then g(X nA T [t M]c ) is a positive supermartingale. Note that 
|a -1 lnxl" -1 = ^(j:) 1 -^ for x < xq and a;" -1 = g(x) 1_ ° for a; > xq. Thus with 
some constant (3 > 

£[<?(X n+1 ) - .g(X„) | X n = a:] < -(3{g{x)) l -± £ [e, M], (52) 

whence 

^ob(^(r.+ l)AT [t , M] c) - 3(^i>AT [e , M ]c)l < -^Xolsftl'^Vlai]^")!' ( 53 ) 

Assume a < 1 so that x— >x l ~ a is decreasing. Then for any sequence p„ > 0, 

EXoigiXn) 1 ^* l {T[tM]c > n} ] > Pn a E Xo [l{g(X n )<p n }l{T [C}M] o>n}} 
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>Pn a {PXo{T[e,M] c ^ n > ) 

so that finally we deduce 

ftir°Pxo(T letM] .>n) (54) 
< E Xo [g(X nAT[CiM] c) - 3(X (n+ i) AT[eM]c )] + (3p n a E Xo {g{X n )l {T[t M]a > n} ). 

Observe that by lt52"|) Ex (g(X n )l^ T M]C >n}) < g(Xo). Then summing the 
inequalities (S^ over n = 0, 1, . . . , m yields 

m m 1 

^X^' 1 " P -^o( T [e,M]= > n) < 3(X ) - E Xo g{X m+1AT[eM]C ) + (3g{X Q )^2p n « 

m 

< 5 (X ) + /3 5 (X ) Vm=l,2,.... 

n=0 

Letting m-KX) we conclude that 

i 

ri>0 n>0 



^2pn a < OO =^ °Px (r[ e ,M]c > 71 ) < 00. 



Now it remains to take p% = n(ln 71) with £ > small enough to get the right 
asymptotics of the tail of Tj e jvf] c since for any non negative random variable T 
and any a > 

E[T a ] < 00 & na ~ lp ( T > 

n) < 00. 

n>0 

If a > 1, we set Y"„ = g(X„)2k so that l|52l) becomes 

E[Yn+i - Yn a I = v] < -PY^ 2 V > L, (55) 

with the special choice of e, M such that L = g(e)^ = g(M)^ (which we can 
always do up to take e smaller or M larger). Hence, by Theorem 1 of [Q, we 
deduce, since {Y n > L} C {T] c > n}, that for a > 1, there exists c < oo 
such that 

£[T£ M] c] < cff(Xo). 



8.4 Proof of Theorem [Ql 

The results <|17fl and (|18|l follow from the estimates on passage-time moments for 
nonnegative stochastic processes of pQ. These are generalizations of Lamperti's 
results for countable Markov chains • 
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Consider the Markov chain on [0, oo[ 

Z n +i = (Z„ - hi/(e Zn ) - F„)l{2 7i _i n y( e -z)_y ?i> o ! z„>o} 

that starts at Zq > and has as an absorption state. Let Sa = inf{n > : 
Z n S A c } the time to escape from the subset A £ [0, oo). We will show that for 
any p < 1/2 and any So > 0, there exists e(6) > such that for any < e < e(S) 

EzoS^^^OiZ 2 ^ 5 ), Z ^oc (56) 
and also that for any p > 1/2 and any e > 

^-]n«,oo[ = °°- (57) 

Then the change of variables Z n = — In X n implies immediately the statements 
1)17(1 . 1(18JI of the theorem as we have T e = S]_i ne)00 [ through this change. It 
follows from (|17|) that T £ < oo a.s. for any Xq < e. This fact entails the 
recurrence of the MC proceeding exactly by the same arguments as in the case 
/(0) > 1. The recurrence is null by l|18|) . The statement l|19|) is a direct 
consequence of ((17)) and ifTBl) . 

To see (|3o]l and (|57|) . we denote by 

f \ _ j E(Z n+ i - Z n \Z n = z) if r = 1, 
^ lzJ ~ \ - Z n |«-|Z„ = z) if r ^ 1. 

By the assumptions EY = 0, E\Y\ 2+Sl < oo we have P(Y > z) = 0(z- 2 - Sl ) 
EY1 Y<Z = 0{z^ 1 - Sl ), EY 2 1 Y<Z = 0{z- Sl ) as z -> oo. By the asumtpion 
/(x) = 1 + o(| lnajl -1-152 ) as a; — » with some 82 > 0, we have f(e~ z ) = 
1 + o(z~ 1 - 52 ) as z -> 00. Then 

Ml( z ) = ^((-^/(e^ 2 ) - ^)l{Y< 2 -ln/(e — )} - Zl{Y>*-la/(e-*)}) 

= O^- 1 - 52 ) + 0(z- 1 -' 51 ), z^oc 

fl 2 (z) = E((-lnf(e~ z ) - Y) 2 l {Y< z-lnf(e-*)} + Z 2 l {Y> z-\n f(e~*)}) 

= 0(z- 2 ~ 2S2 ) + 0(z-' 2 - Sl - S2 ) + EY 2 + 0(z~ 51 ), z-^00, 
and for any 2 < r < 2 + S 2 

fi r .(z) < 2>(1) + \Y\ r ) + z r O(z- 2 - Sl ) = o{z r - 2 ), z~>oo. 

Then for any pel 

2zm(z) + (2p - l)ia{z) = {2 P - l)E{Y 2 ) + 0(z- s ) 

with some 6 > 0. Consequently, Propositions 1 and 2 p. 957 in pQ apply to the 
MC Z n and prove l|56[) and (|57|l concluding the proof of the theorem. 
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8.5 Proof of Theorem 16.11 

Clearly, the transience in cases (1), (2) and (3) follows from (12. '{1 and l|26[) . We 
restrict ourselves to (|23|) . the second case being symmetrical. 
Let us construct the function 

g{X\X 2 ) = (ra InX* - (n - e) lnX 2 )- 1 V cT 1 (58) 

with some fixed d < M. We use the notation X — (X , X 2 ) for shortness. We 
shall prove that there exists (3q > such that for all (3 > (3q 

E( 9 (x n +i) - g(x n ) I g(x n ) = /r 1 ) < o. (59) 

Then for all n > 1 

Eg{X nhT M) < Eg(Xo) = (r a hiX,, 1 - (n - ejln**)- 1 . 

Assume now the contrary of (|23[) : P(t ai < oo) = 1. Then X nAr M — ► X t m a.s. 
as n — ■+ oo. Then by Fatou's lemma 

£<?(X t m) < lim Eg(X nAT M) < Eg(X a ) = {r 2 biX* - (n - ejlnjtf)- 1 . 

n — >oo 

This contradicts the definition of r M as g(X nAr M) > M _1 . 

The proof of (|23|) is reduced now to (|59|l . Let us note that under condition 
r 2 lnX* — (ri — e) In X% = /3 we have: 

r 2 In Xl +1 - (n - e) In X 2 +1 = /? + + h 2 X 2 n + h + Y n (60) 

where 

hi = -r 2 an + (n - e)ai2 > 0, /i 2 = -r 2 ai 2 + (ri - e)a 22 , /i = er 2 > 0, 

r„ = r 2 y n x - (n - e )y„ 2 . 

Our assumptions on and <Zjj for i,j = 1,2 in this case ensure the following 
property: for any 5 > there exists (3q > such that for aJi X 1 ,^ 2 > with 
r 2 lnX 1 - (n -e) lnX 2 = /3 > A) 

/iiX 1 + h 2 X 2 > -5. (61) 

In fact, if condition (1) holds, then also h 2 > 0. Then trivially this quantity is 
non- negative. Otherwise, if h\ > but h 2 < and r 2 < r±, we can rewrite 

^X 1 + h 2 X 2 = h ie lnXl + ^e-ZJ/m+Im/Cn-^-ilin^+fcx* 

If X„ — > oo, the first term with positive coefficient hi > dominates this sum 
due to the fact that r 2 /(ri — e) — 1 < 0. Hence, we have the desired estimate 
for all X 1 > C and all j3 > 0. It remains to chose /3o($) > large enough to 
ensure l|61|l for < X 1 < C. For the rest of the proof we fix S = h/2. Thus 

hiX* + h 2 X 2 + h > h/2 > V/3>/? . (62) 
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The left-hand side of l|59[) equals l\ + /J where 

= E[{d~ l - [3~ 1 )l{Y, l <d-f3-h 1 xi-h 2 x2-h}} 
< d- 1 P(Y n <d-f3-h 1 X^-h 2 X^-h) 

K 

< 



1} 



(/3 + hiXl + h 2 Xl +h- df ln(/3 + h x Xl + hiXl + h-d) 



(63) 



for some k > as E\Y\ 2 < oo. Here we used the fact that (3 is big enough so 
that (3 - d + (h/2) > 0. We have 



r2 

i /3 



E[(([3 ~ h 1 Xj l - h 2 X 2 - h - Y n ) 1 - /3 1 )l{y ri>rf _^_, ll xi -h 2 x2 -h}] 



(/iiXi + /i 2 X2 + h + /3)(/iiXi + fe 2 !2 + /3) 
1 



L {Y" 7l >(i-/3-/iiXi-^2^J-'i} 



03 + hxXl + h 2 Xlf 

Jf3 + h 1 X^ + h 2 X 2 )(h 1 Xl + h 2 X 2 ) 



l{Y n >d-l3-h 1 Xi-h 2 X2-h}] 



(3 

+ h Y X\ + h 2 Xl 
(3 + h x X l n + h 2 Xl + h 

r Y n (c + Y n ) 

X % + h Y X\ + h 2 Xl + h + Y n hYM-H^-h^-H) 

One can bound it from above by 



< 



- {h 1 X} l + h 2 Xl + h + l3)(h 1 X} l + h 2 Xl + 13) 

x ( - h + ftP(y„ < d - 13 - h/2) + +E[Y n l {Yn<d _ _ h/2 



-1 



08 + fell' + ^2^)2 



(64) 



(65) 



- /i/2 
x E 



(3 + hiXl + h 2 X 2 n + h 
Y n (c + Y n ) 



(3 + h x Xl + h 2 Xl +h + Y n 



^{Y n >d-fi-h 1 X^-h 2 X^-h} 



where in the first estimate we used the fact that EY n = 0. Here P(Y n < 
d — (3 - h/2) — > and EY n l{ Yn <d-fi-h/2} ~ * as (3 — > oo since E\Y\ 2 < oo. 

Finally the sequence of the r.v. Z m = „ L+Y l{y„>rf-m} converges to 
zero a.s. as m — > oo and \Z m \ < d~ l \Y n {c + Y n )\ where E\Y n {c + Y n )\ < oo. 
Then by dominated convergence theorem EZ m — > as m — > oo. Since by l|61|l 
13 + hiXl + h 2 Xl + h>/3 + h/2 



E 



Y n (c + Y n ) 



(3 + h x Xl + h 2 Xl + h + Y n 



L {Y n >d-/3-h 1 X} l -h2X^-h,} 



0, (3 — > OO 
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uniformly for all X^,X 2 > with r 2 X^ — (n — e)X^ — (3- Combining these 
facts, we see that the sum of the estimates of (l6*5l) and (l6*3*j) equals ((3 + fti A^ + 
h 2 Xl + h)- 2 (o(l) - ft + o(l) + ft/2 + o(l)) as /3 -> oo uniformly for all X*, X 2 n 
with r2X^ — (fx — z)X\ = /?. Then l|59|l is satisfied and the first part of the 
theorem is proved. To prove (|24[1 . we notice that 

Xl +l < Xle-^^+Yl (66) 

Xl +l > X 1 n e- a ^- d( - M ^ x ^+ r ^ onr M >n (67) 

with d(M) = ai2e~ ^i-' M and S — . We therefore can prove that T^X 1 ) 
satisfies the analogue of the bound (|28fl of Theorem 12.11 by (|66fl (note that we 
only used in the proof of (|28|) the upper bound), whereas we can use (|67fl to 
get bounds Ijllfl and (|12|l on T^X 1 ) under condition t m = oo. In fact, we can 
apply the same arguments than in the proofs of Theorems 12.11 14.11 and 14.21 by 
considering the [X n , X nATM )„> instead of (X n ,X n ) n > yielding Ex^T^X 1 ) A 
Tm] < oo and £ , x [2l(^ 1 )] < °o (note here that in the neighborhood of the 
origin, the correction d(M)(X^) s is small). Then as in one-dimensional case, 
we can derive that Px (T[ e ^] C (X 1 ) A tm < oo) = 1< 

Acknowledgment We are very grateful to O. Zeitouni for pointing out some 
useful references. 
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